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We use our recently proposed microscopic modeling in the framework of linear response theory to 
investigate the complete phonon dispersion, the phonon density of states, certain phonon-induced 
electronic charge distributions and charge fluctuations (CF's) for anomalous soft modes of metal- 
lic BaBiC>3 in its simple cubic phase where superconductivity with T c up to 32 K appears. The 
theoretical approach already has been applied successfully to the cuprate high-temperature super- 
conductors (HTSC's), simple ionic crystals (NaCl, MgO) and perovskite oxides (SrTiC>3, BaTiOa). 
It is well suited for materials with a strong component of ionic binding and especially for "ionic" 
metals. In particular, the giant phonon anomalies related to the breathing vibration of the oxygen 
as found experimentally in superconducting doped Bao.6Ko.4Bi03, resembling those observed in 
the high T c cuprates, are investigated. The origin of these anomalies is explored and attributed 
to a strong nonlocal coupling of the displaced oxygen ions to CF's of ionic type, essentially of the 
Bi6s- and Bi6p orbital. This points to the importance of both of these states at the Fermi energy. 
Starting from an ab-initio rigid ion model (RIM) we calculate the effect on the lattice dynamics and 
charge response of the most important electronic polarization processes in the material, i. e. CF's 
and dipole fluctuations (DF's). Taking into account these electronic degrees of freedom in linear 
response theory, we obtain a good agreement with the measured phonon dispersion and in particular 
with the strong phonon anomalies. 
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I. INTRODUCTION 

Potassium doped Bai^^K^BiOa with a transition 
temperature T c up to 32 K for x rs 0.4 has been in- 
tensively studied with no final resolution of its high- 
temperature superconductivity. Different from the 
cuprate high-temperature superconductors (HTSC's) 
its bonding is three dimensional and no antiferromag- 
netic ordering exists for the insulating parent com- 
pound indicating that there are no strong electronic 
correlations near x — 0. Common to both ma- 
terials, however, is the strong component of ionic 
binding favoring long-ranged Coulomb interactions 
and a related strong nonlocal electron-phonon cou- 
pling in terms of localized charge fluctuations (CF's) 
of the ionic orbitals. This type of coupling has 
been shown in the cuprate-HTSC's to lead to the 
generic phonon anomalies of the high-frequency oxy- 
gen bond-stretching modes (OBSM)^^^. Recently, 
giant phonon anomalies for similar OBSM (oxygen 
breathing modes) also have been found by inelastic 
neutron scattering measurements in superconducting 
Bao.6Ko.4Bi03 7,8,9 . Thus, while spin-degrees of free- 
dom seem to play no decisive role for superconduc- 
tivity in Ba-Bi-0 the effects of nonlocal coupling of 
lattice and charge degrees of freedom in terms of 



ionic CF's are still present, likewise as in the cuprate- 
HTSC's, and provide a coupling channel for pairing in 
a phonon-mediated mechanism. 

Bai-^K^BiOa has a complex structural phase 
diagram^. For the insulating composition < x < 
0.12 the structure is monoclinic and may be de- 
rived from the cubic perovskite structure by simul- 
taneous octahedral tilting (t) and symmetric oxygen 
breathing-mode (b) distortions (frozen-in OBSM in- 
stability at the i?-point of the cubic phase). For 
0.12 < x < 0.37 the structure is still insulating but 
orthorhombic with a t, but no b distortion and for the 
superconducting, metallic composition (0.37 < x < 
0.53) the structure is simple cubic, i. e. no t and b 
distortion. 

It is widely believed that the insulating phases are gen- 
erated by charge density instabilities associated with 
breathing and tilting distortions but it has proven dif- 
ficult to establish this by first principles calculations. 
Recent investigations^ indicate, in contrast to earlier 
results^, that the local-density approximation (LDA) 
seriously underestimates the breathing distortions and 
yields a metallic ground state with no sign of a charge- 
density-wave instability. 

Moreover, first principles linear response calculations 
of the phonon dispersion of cubic Ba-Bi-O^i overes- 



2 



timate the high-frequency branches by about 20 % in 
comparison with the experiments. In particular the 
most strongly renormalized anomalous OBSM, Of, 
at the i?-point, which according to the experiments^ 
is found at about 10.5 THz has a frequency of 15.7 
THz in the calculations presented in Ref. [ll|. For the 
modes with lower frequencies there is a better agree- 
ment of the order of 10 %. Thus, linear response calcu- 
lations based on LDA-like electronic band structures, 
using a virtual-crystal approximation seem to under- 
estimate the kinetic part of the charge response, i. e. 
the electronic polarizability matrix. In order to inves- 
tigate this point, we study different cases for the po- 
larizability matrix H KK > (q ) in our microscopic model 
approach to extract the most important electronic de- 
grees of freedom for the charge response. The indices 
k, k! denote the orbital degrees of freedom in an ele- 
mentary cell of the crystal. 



Empirical shell model calculation; 
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incorporating 



metallic screening by the Lindhard function are also 
unable to reproduce the downward dispersion and 
softening of the phonon anomalies, instead an upward 
dispersion is calculated. Finally, calculations in Ref. 
13l based on a tight-binding representation of the elec- 
tronic contribution to the dynamical matrix— and an 
empirical contribution for the short-ranged force con- 
stants fitted to the experimental frequencies lead to 
a remarkable softening for several modes especially 
around the M- and the i?-point with strongest soften- 
ing for the (1, 1, 0) direction, in disagreement with the 
experiments where only the OBSM of breathing type 
are strongly renormalized and the largest softening is 
along the (1, 1, 1) direction. 

In the present work we calculate within our micro- 
scopic modeling of the electronic charge response in 
linear response theory the complete phonon disper- 
sion, the phonon density of states, the phonon-induced 
electronic charge redistribution and the CF's for the 
OBSM, Of, Of, Of, of simple cubic Ba-Bi-O. In the 
past, this theoretical approach, which is well suited 
for systems with a strong component of ionic bind- 
ing, has been applied successfully to the cuprate- 
HTSC'si^££, i.e., "ionic" metals in which the 
band crossing the Fermi level has admixture of both 
anion and cation orbitals, simple ionic crystals^ and 
the perovskite oxides (SrTi0 3 , BaTi0 3 )^. 
From a general point of view our treatment of the elec- 
tronic density response and lattice dynamics in terms 
of dipole fluctuations (DF's) and CF's can be con- 
sidered as a microscopic (semi-ab initio) implemen- 
tation of the phenomenological dipole-shell model or 
the charge-fluctuation models, respectively. For a gen- 
eral formulation of phenomenological models for lat- 
tice dynamics that use localized electronic variables 



as adiabatic degrees of freedom, see, for example, Ref. 
171 . This formulation covers shell models, bond-charge 
models and charge-fluctuation models. While in such 
an approach the coupling coefficients are treated sim- 
ply as empirical fitting parameters the essential point 
in our scheme is that all the couplings are microscop- 
ically well defined and can be calculated. 
The ionic nature of Ba-Bi-0 is described by an 
ab initio rigid ion model (RIM) leading to a local 
rigid charge response and electron-phonon interaction 
(EPI), respectively. The missing nonlocal, nonrigid 
part of the electronic density response and EPI is 
expressed by microscopically well defined CF's and 
DF's on the outer shells of the ions. Starting from 
the RIM as a reference system the effect of the non- 
rigid electronic polarization effects in terms of CF's 
and DF's on the phonon dispersion and, particularly, 
the anomalies is calculated and compared with the 
experiment. 

The article is organized as follows. In section 2 we 
outline the theory and modeling to provide a better 
reading of the article. Section 3 presents our calcu- 
lated results of the phonon dispersion and gives a de- 
tailed discussion of the effects of screening by DF's 
and CF's on the phonon anomalies by comparing the 
corresponding results with those as obtained from the 
RIM as a reference system. Moreover, phonon density 
of states, the charge redistribution and the CF's for 
the OBSM of breathing type are presented. Finally, 
section 4 contains a summary and discussion. 



II. SKETCH OF THE THEORY AND 
MODELING 

In the following a survey of the theory and modeling is 
presented. A more detailed description can be found 
in Ref.Q and in particular in Ref. 15J where the calcu- 
lation of the microscopic coupling parameters is given. 
The local, rigid part of the electronic charge response 
and the EPI is approximated by an ab initio RIM 
taking into account ion-softening in terms of (static) 
effective ionic charges and scaling of the short-ranged 
part of certain pair potentials between the ions to sim- 
ulate covalence effects in the model. This is done in 
such a way that the energy-minimized structure is as 
close as possible to the experimental one-^. Structure 
optimization and energy minimization is very impor- 
tant for a reliable calculation of the phonon dynamics 
through the dynamical matrix. 

The RIM with the corrections just mentioned serves as 
an unbiased reference system for the investigation of 
the effect of the nonrigid electronic polarization pro- 
cesses on the phonon dynamics. The latter are mod- 
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eled in form of electronic CF's on the outer shells of 
the ions. Especially in the metallic state the CF's 
dominate the nonlocal contribution of the charge re- 
sponse and the EPI. In addition DF's are admitted 
in our approach£*!£. Thus, the basic variable of our 
model is the ionic density which is given in the per- 
turbed state by 

Pa{r, Qx,p a ) = P° a (r) + Q\Px ¥ {r) + p a • rpg(r). 

A 

(1) 

is the density of the unperturbed ion, as used in the 
RIM, localized at the sublattice a of the crystal and 
moving rigidly under displacement. The Q\ and p^ F 
describe the amplitude and the form-factors of the 
CF's and the last term in Eq. (p]) represents the dipo- 
lar deformation of an ion a with amplitude (dipolc 
moment) p Q and a radial density distribution p®. r 
denotes the unit vector in the direction of r. The p^ F 
are approximated by a spherical average of the orbital 
densities of the ionic shells calculated in LDA taken 
self-interaction effects (SIC) into account. The dipole 
density p„ is obtained from a modified Sternheimer 
method in the framework of LDA-SIC 15 . All SIC cal- 
culations are performed for the average spherical shell 
in the orbital-averaged form according to Ref. LJjJ For 
the correlation part of the en ergy per electron, e, the 
parametrization given in Ref. [19| has been used. 
The total energy of the crystal is obtained by assum- 
ing that the density can be approximated by a super- 
position of overlapping densities p a . The p° a in Eq. 
([1]) are also calculated within LDA-SIC taking envi- 
ronment effects, via a Watson-sphere potential and 
the static effective charges of the ions into account. 
The Watson-sphere method is only used for the oxy- 
gen ions and the depth of the Watson-sphere potential 
is set as the Madelung potential at the corresponding 
site. Finally, applying the pair-potential approxima- 
tion we get the total energy 



E{R,Q=Y, E ^0 + \ 

a.a ~ 



E 

(a,a)#(b,/3) 



$ Q/3 (r£-r»,c). 

(2) 

The energy E depends on both, the configuration of 
the ions {R} and the electronic (charge) degrees of 
freedom (EDF) {£} of the charge density i.e. {Q\} 
and {p Q } in Eq. |T]). E* are the energies of the single 
ions, a, b denote the elementary cells and a, (3 the 
corresponding sublattices. The second term in Eq. 
@ is the interaction energy of the system expressed 
in terms of anisotropic pair interactions <S> a p- Both 
E* and <& a p in general depend upon £ via p a in Eq. 
©■ 



The pair potentials in Eq. $2$ can be separated into 
long-ranged Coulomb contributions and short-ranged 
terms as follows: 



*a/3(R, C) 



ZryZ 



R 



{Z a pp + Z[3p a ) ■ — -L 



R 3 



3 ( Pa .R)(R.p,) + ^ (R;a (3) 



$ a/3 (R, C) =K a U p {R, C) + KpU a (IL, c) 

+ W^(R,C) + G^(R,C)- (4) 

The first term in Eq. ([3]) describes the long-ranged 
ion-ion, the second the dipole-ion and the third and 
fourth term the dipole-dipole interaction. Z a and Zp 
are the variable charges of the ions in case CF's are 
excited. The latter reduce to the ionic charges for 
rigid ions. K a and Kp are the charges of the ion 
cores. $ Q( 3 represents the short-ranged interactions. 
These can be expressed by the following integrals 



U a (TL,() 



d 3 rp a (r, C) 



Rl 



l 

R 



R 



R 3 



(5) 



W af) (R,0 = I d 3 rl dV[p a (r,C)Mr',C)x 

1 1 (r + r') • R 

|r - r' - R| i? R 3 

G Q/3 (R,C) = J d 3 r[p af3 (r,Oe(p af3 (T,0) 

-p Q (r,C)e(p tt (r,C)) 
-Mr-RCMMr-RC))], 



(G) 



(7) 



with 



P*p(r, = p a (r, C) + Pp(r - R, 0- 



(8) 



K a Up (R, C) yields the short-ranged contribution of 
the interaction between the core a and the density 
pp according to Eq. (JTJ) . W^^R, () represents the 
short-ranged Coulomb contribution of the interaction 
of the density p a with the density pp and G a p (R, Q) is 
the sum of the kinetic one-particle- and the exchange- 
correlation (XC) contribution of the interaction be- 
tween the two ionsi^. The short-ranged part of the 
potentials and the various coupling coefficients are cal- 
culated numerically for a set of distances R between 
the ions. The results so obtained are then described 
by an analytical function of the form 



f(R) = ±exp(a + /3i? 



(9) 
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a, j3, and 7 in Eq. © are fit parameters. 
From the adiabatic condition 



dE(R,Q 
8( 







(10) 



an expression for the atomic force constants, and ac- 
cordingly the dynamical matrix in harmonic approxi- 
mation consistent with linear response theory can be 
derived: 



t,f(q)= C(q) 



pet 

'ij 



RIM 



1 



y/M a M k r , 

(11) 

The first term on the right hand side denotes the con- 
tribution from the RIM. M a , Ms are the masses of 
the ions and q is a wave-vector from the first Bril- 
louin zone. The quantities B(q) and C(q) in Eq. 
(fTT|) represent the Fourier transforms of the electronic 
coupling coefficients as calculated from the energy in 
Eq. ([2|), and the pair potentials in Eqs. (J3|)-([8|), re- 
spectively: 



II -1 contains the kinetic single particle contribution 
to the interaction C and V the Hartree and exchange- 
correlation contribution. C _1 needed for the dynam- 
ical matrix and the EPI is closely related to the (lin- 
ear) density response function (matrix) and to the 
inverse dielectric function (matrix) e -1 , respectively. 
Comparing with calculations of the phonon dispersion 
and the EPI using the linear response method in form 
of density functional perturbation theory within local 
density approximation (LDA), these calculations cor- 
respond to calculating II and V in DFT-LDA. On the 
oder hand, in our microscopic modeling DFT-LDA- 
SIC calculations are performed for the various den- 
sities in Eq. |T]) in order to obtain the coupling co- 
efficients B and V. SIC as a correction for a single 
particle term is important for contracting in particu- 
lar localized orbitals. Written in matrix notation we 
get for the density response matrix the relation 



n(i + yn)- 



e = 1 + VTI. (16) 



The CF-CF submatrix of the matrix II can for exam- 
ple approximated from a tight-binding approximation 
of a single particle electronic bandstructure. In this 
case the electronic polarizability II reads 



Bab 
KB 
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d 2 E(R,Q 



(12) 
(13) 



k denotes the EDF (CF and DF in the present model) 
in an elementary cell. The B coefficients describe 
the coupling between the EDF and the displaced ions 
(bare electron-phonon coupling), and the coefficients 
C determine the interaction between the EDF. The 
phonon frequencies, £J ff (q), and the corresponding 
eigenvectors, e a (qtr), of the modes (q<r) are obtained 
from the secular equation for the dynamical matrix in 
Eq. (IT>, i.e. 



E 

3,3 



^(q)ef(q) 



^ 2 (qK(q) 



(14) 



IW(q, uj = 0) 



- V 

N 



/n'(k + q)-/ n (k) 
£„/(k + q) -£„(k) 



X [C*„(k)C Kn ,(k + q)] [C:,„(k)^„-(k + q)]* . 

(17) 

/, E and C in Eq. (fl~7|) are the occupation numbers, 
the single-particle energies and the expansion coeffi- 
cients of the Bloch-functions in terms of tight-binding 
functions. 

As a measure of the strength of the EPI the self- 
consistent change of an EDF on an ion, induced by 
a phonon mode (q, a), can be derived in the form 



-^X KQ (q)u a ( q( T) 



e ^ R 5 ee <5C K (qa)e^ Ra 
(18) 



with the displacement of the ions 



The lengthy details of the calculation of B and C can- 
not be reviewed in this paper. They are given in Ref. 
15. In this context we remark that the coupling matrix 
C KK > (q ) of the EDF-EDF interaction, whose inverse 
appears in Eq. (fTTj) for the dynamical matrix, can be 
written in matrix notation as 



<( q( x) 



2M a uj a (q) 



1/2 



e a (qa)e iqRa = u Q (qa)e lclRa 
(19) 



The quantity X in Eq. fTS]) . i. e. the self-consistent re- 
sponse per unit displacement of the EDF is calculated 
in linear response theory as 



c = ir 



v. 



(15) 



X(q) = n(q) £ - 1 (q)B(q) = C7- 1 (q)B(q). (20) 
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Another measure of the EPI for a certain phonon 
mode (qer) is provided by the change of the self- 
consistent potential in the crystal felt by an electron 
at a space point r, i.e. 6V c g(r, qcr). Averaging this 
quantity with the corresponding density form factor 
p K (r — R*) of the EDF located at R*, we obtain 

6V*(qa) = J dVp K (r - B.%)8V eS (r, qcr). (21) 

This gives an orbital resolved measure for the strength 
of the EPI in the mode (qcr) mediated by the EDF 
considered. For an expression of 6V* (qcr) in terms of 
coupling coefficients in Eqs. (fl~2|) and (Tl3|) . see Refi 2 ^. 

III. RESULTS AND DISCUSSION 
A. Results within the reference System (RIM) 

The rigid, local contribution of the electronic charge 
response and EPI is approximated by an ab initio RIM 
with corrections for covalence effects in terms of effec- 
tive ionic charges (ion softening) and scaling of the 
short-ranged pair potential between the Bi and oxy- 
gen ion. In this way we obtain a suitable reference 
system, see Fig. Q]for the calculated dispersion, which 
subsequently allows for the investigation of the charac- 
teristic nonrigid screening effects in terms of CF's and 
DF's on the phonon dispersion and charge response. 
The covalent corrections are important because a cal- 
culation of the phonon dispersion of cubic Ba-Bi-0 
with nominal charges leads to a significant overesti- 
mation of the width of the phonon spectrum. More- 
over, various unstable modes down to -12 THz emerge 
if imaginary frequencies of unstable modes are rep- 
resented as negative numbers. The mode with the 
highest frequency (« 35 THz) in such a model is the 
oxygen- breathing mode, Og, at the R point which, 
on the other hand, according to recent experiments 9 
is anomalously soft (~ 10.5 THz), see Fig. [T] Similar 
results are found in calculations of the phonon dis- 
persion of cuprate high-temperature superconductors 
when nominal ionic changes are usedi^. In our calcu- 
lation we find a set of effective (static) ionic charges 
(Ba 1.5+, Bi 2.7+, O 1.4-) and a covalent scaling 
of the Bi-0 pair potential which leads to an energy 
minimized structure in full agreement with the ex- 
perimental one for cubic Ba-Bi-0 (a — 4.2742 A)^. 
As a general rule, partial covalence reduces the ampli- 
tude of the static ionic charges in mixed ionic-covalent 
compounds like Ba-Bi-O, because the charge transfer 
from the cations to the anions is not complete as in 
the entirely ionic case. 



The phonon dispersion based on this set of reduced 
ionic charges is displayed in Fig. [T] (model 1). Potas- 
sium doping is taken into account using a virtual-mass 
approximation and implicitly via the reduced static 
charges. The phonon dispersion calculated with such 
a reference model means a large improvement com- 
pared to the model with nominal charges (not shown). 
In particular, the width of the spectrum is consider- 
ably reduced and compares well with the experiments 
if the missing renormalization related to the nonrigid 
charge response (CF's and DF's) is taken into account, 
see subsection MI BI Moreover, unstable modes only 
appear at the R and the M point. The mode at R 
is triply degenerate and represents rotations around 
any cubic axis. The rotation is opposite in neighbour- 
ing elementary cells in all cubic directions. The other 
unstable mode at M is similar to the most unstable 
R point mode, with the exception that the rotation 
of the octahedra is in the same sense in adjacent cells 
along the z-axis. Note in this context that, as men- 
tioned in the introduction, the structural phase tran- 
sitions in K-doped Ba-Bi-0 are characterized by rota- 
tions (tilts) of the BiOe octahedra around several dis- 
tinct axisi^ and in cubic Bao.6Ko.4Bi03 the frequency 
of the rotational R point mode is extremely low (Fig. 
[T]), consistent with the tendency of soft rotational 
modes at R in harmonic theory. Thus, our results 
within the harmonic approximation highlight the im- 
portance of the ionic forces for the tilt instabilities and 
indicate that anharmonic contributions should stabi- 
lize these modes. Frozen-phonon calculations support 
this argument because they predict highly anharmonic 
potential behaviour for the tilting of the BiOg octa- 
hedra corresponding to the fi-poin t 11 ' 22 . 
In SrTi03 the rotational modes at R also are found 
most unstable in our calculations^ corresponding to 
an antiferrodistortive transition in agreement with 
the experiment while our calculation for BaTi03 find 
the ferroelectric mode at V as the most unstable 
one again in accordance with the experimental sit- 
uation. Interestingly, we also obtain in case of the 
HTSC La2CuQ 4 18 i 20 one partially unstable branch 
with the tilt mode at the A-point. Frcezing-in of 
this distortion points correctly to the experimentally 
observed structural phase transition from the high- 
temperature tetragonal to the low-temperature or- 
thorhombic structure. Again this rotational instabil- 
ity is already present in the RIM and thus also driven 
by the long-ranged ionic interactions likewise as in Ba- 
Bi-0 and Sr-Ti-O. In contrast to the situation in Ba- 
Ti-0 the strongly polar ferroelectric TO mode at T, 
where the oxygen anions vibrate coherently in the op- 
posite direction to the Bi ion, is stable in Ba-Bi-O. We 
get 8.6 THz in the RIM with effective ionic charges 
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FIG. 1: Calculated phonon dispersion of Ba-Bi-0 in the cubic perovskite structure in the main symmetry directions 
A ~ (0, 0, 1), E ~ (1, 1, 0) and A ~ (1, 1, 1) for the ionic reference model with covalent corrections (RIM, model 1). 
The open symbols represent the experimental results for cubic Bao.6 K0.4 BiO^. The full squares are new experimental 
results of the same material. The arrangement of the panels from left to right according to different irreducible 
representations is as follows: |Ai (— ), A 2 (• • • )| A 3 (— )|Ei (— )|E 2 (— ), S3 (• - - )| |S 4 (— )|Ai (— ), A 2 (■ ■ ■ )|A 3 (— )|. The 
measured anomalous soft oxygen-breathing vibrations at the X, M and R point show the following frequencies: Og ~ 
12.8 THz, Og ~ 11.8 THz and Og ~ 10.5 THz, respectively. Imaginary frequencies of unstable modes are represented as 
negative numbers. 



and 9.7 THz in the RIM with nominal ones. Not only 
the frequency of the ferroelectric mode is significantly 
reduced when passing from the model with nominal 
ionic charges to the model with effective charges but 
also the LO-TO mode splittings. Inspection of Fig. 
[T] with respect to the oxygen-breathing modes, dis- 
played in Fig. [21 we find that instead of being soft 
as in the experiment the latter have the highest fre- 
quency at the X (Of), M (0%) and R (0§) point, re- 
spectively. Consequently, the characteristic downward 
dispersion of the anomalous phonon-branches cannot 
be described at all within the RIM, i.e., by a rigid, 
local charge response. In the following we will show 
that the anomalous dispersion and phonon softening 
characteristic for cubic, metallic Ba-Bi-0 is essentially 
related to strong nonlocal, nonrigid EPI effects medi- 
ated by CF's of the Bi6s and Bi6p orbitals, pointing 
to a substantial character also of Bi6p electrons in the 
electronic state near the Fermi energy. 
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FIG. 2: Displacement patterns of the anomalous oxygen- 
breathing modes Og, Og 1 and Og from left to right for 
model 3. Frequencies are given in units of THz. 



B. The effects of a nonrigid, nonlocal charge 
response via CF's and DF's 

Starting from the reference systems (RIM) of the last 
subsection DF's and CF's as nonrigid electronic de- 
grees of freedom arc additionally allowed for the mod- 
eling. This means that on each ion the electrons can 
redistribute under atomic displacements in such a way 
that dipole and charge fluctuations are induced on 
that ion in order to minimize the energy. 
In Fig. [3] we display our calculated results of the 
phonon dispersion taking additionally DF's into ac- 
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FIG. 3: Same as in Fig. Q]with the calculated results from the model described in the text where additionally to the 
RIM (model 1) dipole fluctuations are taken into account (model 2). 



count (model 2). In this calculation the ab initio val- 
ues of the dipole polarizability a as calculated from 
the Sternheimer method in the framework of DFT- 
LDA-SIGi^ for the corresponding single ions are re- 
duced by 50 % which leads to a better agreement of 
the phonon dispersion with the experiment. Taking 
the ab initio values for a as calculated for the single 
ions with the Sternheimer method overestimates the 
dipole polarization in the crystalline environment and 
leads to TO-frequencies being to low as compared to 
the experiment. These findings also hold true for our 
calculations of dipole screening in the cuprate based 
HTSC's^^ where in addition the dipole polarizability 
is found to be very anisotropic dominating along the 
(ionic) c-direction. 

Comparing Fig. [3] with Fig. Q] we find that the 
width of spectrum improves considerably and no addi- 
tional unstable branches emerge. Moreover, the LO- 
TO splittings are further reduced and the TO-modes 
are well described. In particular, the ferroelectronic 
TO-mode at T is reduced from 8.6 to 6.7 THz and 
is now in quite a good agreement with experiment^. 
However, the anomalous branches with the oxygen- 
breathing modes at X, M and R have an increas- 
ing dispersion when propagating from the T point to- 
wards the zone boundary and show no sign of soften- 
ing. Likewise as in the RIM, these modes have the 
highest frequencies in the spectrum. 
Next, we discuss the influence of the ionic CF's on the 
phonon dispersion. The results of a calculation where 



additionally to DF's (model 2) CF's are included are 
shown in Fig. Q] (model 3). Comparing the calculated 
results with the experiments now a good agreement 
is achieved. This is particularly true for the phonon 
anomalies which cannot be described well by other ap- 
proaches up to now as stated in the introduction. First 
principle calculations^ within DFT-LDA for Ba-Bi-0 
taking K-doping into account by a virtual-crystal and 
virtual-mass approximation report that the breathing 
distortions are seriously underestimated in the LDA. 
These calculations yield for the Og mode a frequency 
of 15.7 THz compared to 10.5 THz in the experiment^. 
We argue in the following that a reason for the high 
frequency in this calculation is the electronic polariz- 
ability n KK /(q) from Eq. (fTT)) , i.e. the kinetic single 
particle part of the charge response, which seems un 
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derestimated in LDA-type calculations as in Ref. 
or Ref. [2J, respectively. 

In order to investigate this question in our approach, 
Tl KK > (q ) is modeled in an appropriate way to find out 
which electronic degrees of freedom are most impor- 
tant for the phonon anomalies in Ba-Bi-O. So, we ex- 
press H KK i (q ) in Fourier-transformed (parametrized) 
form 



'(q)=E n °> 4q(R ^ RK) - 



(22) 



In the metallic phase the electronic partial density of 
states (PDOS) at the Fermi level Z k (ef) is related 
to the polarizability matrix in the long wave length- 
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limrfei*& according to 

]TlW(q^O) = Z K ( £F ), (23) 

re' 

and the total density of states at energy e is given by 



(24) 



In our calculations for the metallic phase of the 
cuprate HTSC's we have found^^ that neglecting the 
q-dependence in the polarizability and taking into ac- 
count the diagonal on-site parameters II KK alone al- 
ready provides a good approximation for the calcu- 
lation of the phonon dispersion, in particular for the 
anomalous generic OBSM in these materials. This 
means that the parameters IT K K are approximated by 
the corresponding PDOS. However, such a diagonal 
approximation does not hold for the insulating state 
of a material. Due to the charge gap in the spectrum 
off-site elements 11°* necessarily must occur and in- 
terfere in such a way that the following compressibil- 
ity sum rules for the kinetic single particle part of the 
charge response are satisfiedi^. 



]Trw(q^o) = (%), 



£lW(q-0) =0(q 2 ). 



In Fig. Owe display the calculated effect of the Bi6s-, 
Bi6p- and 02p charge fluctuations for the softening of 
the frequency of the anomalous OBSM Og, Og and 
Og, respectively by varying the diagonal on-site pa- 
rameters IT(Bi6s) = ZBi6s(£F), n (Bi6p) = Zbi6p(£f) 
and n(02p) = Z 02v {e F ). We find that the CF's in 
the Bi6s- and Bi6p orbitals both are important for 
the phonon softening, while the dependence on the 
02p CF's is very weak. The values actually chosen 
for the calculation of the phonon dispersion in model 
3 are n(Bi6s) = l.OOeV -1 , II(Bi6p) = 0.5eV _1 , 
II(02p) = 0.2 eV -1 . The latter value, to a large ex- 
tent irrelevant for the softening, was taken to optimize 
the phonon frequencies for modes where exclusively 
02p CF's are excited. The remarkable softening with 
increasing Zbi6p(£f) indicates that the band crossing 
the Fermi energy consists not only of the Bi6s and 
02p electrons but also has a substantial admixture of 
Bi6p electrons. This is also concluded from electronic 
band structure calculations, see Ref. 11 and refer- 



see Ref. 

ences therein, and is consistent with the static effec- 
tive charge for Bi in our model pointing out a partial 
occupation of the Bi6p orbital. The latter is the most 
delocalized one of the the orbitals used in our calcu- 
lation (see Fig. [6]) and consequently very effective for 
screening and softening of the oxygen bond-stretching 

(25) vibrations. 

Such an enhancement of softening, related to the de- 
localized components of the electronic state at the 

(26) Fermi energy, also has been obtained in our calcula- 



9 












• 


• 1 


V ° 


• o-. 


































▼ ... 














} 1 


2 


3 


4 


5 E 



(a)n(Bi6s) 



(b)n(Bi6 P ) 




(c)n(02 P ) 



FIG. 5: Calculated effect of the Bi6s-, Bi6p and 02p charge fluctuations on the softening of the frequency of the anomalous 
oxygen bond-stretching (breathing) modes Ob, Ob and Ob, respectively, for varying polarizabilities Il(Bi6s), Il(Bi6p) 
and II(02p). II is given in eV _1 . X-point (full line), M-point (dotted line), J?-point (broken line). 
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FIG. 6: Charge density form factors (times Airr 2 ) for the 
Bi 2 ' 7+ and O 14- ion. r radial distance in units of as ■ Bi6s 
(full curve); Bi6p (dotted curve); 02p (broken curve). 



tions of the generic anomalous OBSM in the cuprate 
HTSC's. In this case we have shown that with in- 
creased hole doping the inclusion of the delocalized 
Cu4s state into the orbital basis is an important fac- 
tor for the strength of the anomalies in the optimally 
and overdoped state2&2&. Thus, doping (similar like 
pressure) causes a progressive modification of the or- 



bital character of the electronic state in the cuprate- 
HTSC's. Quite generally, our investigations in the 
cuprates point to an interrelation between electronic 
structure changes at ef, selective phonon softening 
via strong nonlocal EPI of CF' type and, finally, to 
the appearance of phonon-mediated superconductiv- 
ity in an "ionic" metal. A corresponding connection 
seems also to hold in doped Ba-Bi-O. 
A change of the electronic state in Ba-Bi-O, i. e. a 
partial shift of the spectral weight from Bi6s to Bi6p, 
as suggested from our calculations of the OBSM, ob- 
viously occurs parallel with K-doping because the 
OBSM anomalies are not present in the weakly doped 
material. It has been argued in Ref. that K- 
doping creates localized holes on the 02p orbitals, 
rather than on Bi atoms. One might speculate, that 
an increased population of the Bi6p states as the most 
delocalized ones lowers the energy by an enhanced hy- 
bridization with the 02p hole states making the latter 
itinerant at a certain doping level in the metallic state. 
An increasing Bi6p character of the electronic state at 
e f also should lead from a more rounded Fermi surface 
to a more distorted anisotropic shape. Such a change 
of the geometry of the Fermi surface would of course 



be mapped to the "phase space factor" (nesting func- 
tion) and also to the matrix elements of the electron- 
phonon spectral function a 2 F(uj) and, finally, to the 
superconducting T c in a phonon-mediated mechanism. 
Altogether from our results a multi-orbital model in- 
volving CF's, DF's and the corresponding interactions 
V in the charge response, Eqs. (fTSj) . (fTo) , between 
Bi6s-, Bi6p- and 02p states is needed to obtain a reli- 
able description of the phonon dispersion in an " ionic" 
metal like Ba-Bi-O, as shown in Fig. 0J 
Comparing the results of the dispersion in Fig. [3]and|3] 
we recognize that the LO-TO splittings at T are closed 
by the metallic screening due to the CF's. Such a 
screening of the longwavelength polar modes does not 
imply charge localization on a microscopic scale and 
can also be achieved in a homogeneous electron gas 
model. However, localization of charge on a micro- 
scopic scale is necessary for the screening of the zone- 
boundary modes like Of, Og and Og. In this case 
nonlocally induced ionic Bi6s- and Bi6p CF's lead to 
the strong softening. Compare with Fig. [7J where the 
nonlocal part of the phonon-induced charge density 
redistribution is shown for Og, Og and Og. The re- 
lated strong softening of the modes by the nonlocally 
excited nonrigid charge response via CF's and DF's as 
compared to the RIM can be extracted from Fig. [8l 
In context with the vanishing of the LO-TO splitting 
characteristic phonon branches along A, S and A with 
the (stable) ferroelectric TO mode at T (6.744 THz) 
appear in the phonon dispersion which are well de- 
scribed in our model. Even details, like the steep 
slopes near T, resulting from the depression of the 
LO frequency by metallic screening, and the flatten- 
ing out of the branch in A direction is obtained in our 
calculation. 

In order to investigate the question how the phonon 
dispersion looks like, if we use an electronic density 
of states at ep for Ba-Bi-0 which is typical for a 
LDA-like calculation 24 in our model for the polariz- 
ability, we have performed the calculations shown in 
Fig. [S The total density of states from Ref. is 0.8 
e y-i taking for n(02p) the same value as in model 
3 and using the same ratio for II(Bi6s) and n(Bi6p) 
we obtain: n(Bi6s) = 0.4 eV _1 , II(Bi6p) = 0.2eV~\ 
II(02p) = 0.2 eV -1 (model 4). Comparing the results 
for model 4 in Fig. [5] with those in Fig. 0] we find that 
the softening of the oxygen-breathing mode anomalies 
is only poorly developed, consistent with the ab ini- 
tio calculations based on DFT-LDA in Ref. EH. This 
indicates that the ab initio LDA calculations under- 
estimate the EPI for the OBSM. 
The underestimation of the EPI in the breathing 
modes can also be deduced from Table Q] in terms of 
the CF's. Here the CF's according to Eq. (TTg]) of the 
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5C,Bi6s <5CBi6p <5C02p 


of 


39.369 27.534 -4.169 
32.333 18.285 -1.764 


of 


61.329 42.740 -5.791 
46.309 26.227 -2.092 


of 


83.657 57.859 / 
56.807 32.220 / 



Table I: Calculated phonon-induced charge fluctuation 
S( K (qa) for model 3 (first row) and model 4 (second row) 
for Of , Of and Of according to Eq. JT8]). Units are 
in 10~ 3 electrons. 8^ K (qa) > means that electrons are 
pushed away. / means, that no CF's are excited. 

Bi6s-, Bi6p- and 02p orbitals generated by Of, Og 
and Og are listed. We recognize that in our modeling 
of the electronic state with enhanced PDOS for Bi6s 
and Bi6p at ep as compared to the LDA-based result 
the CF's are significantly larger leading to the correct 
description of the anomalies in the phonon dispersion. 
The strong nonlocal EPI of ionic CF type in the 
oxygen-breathing modes is also illustrated in the con- 
tour plots of the nonlocal part of the phonon-induced 
charge density redistribution Sp(r, qcr) for the OBSM 
in Fig. [3 i. e. 

5p(r, qtr) = ]T SQ(qa)p K (r - R*), (27) 

with the CF's SQ according to Eq. (TTg)) and the 
form factors p K for the CF's as in Eq. ([1]). The total 
phonon-induced charge redistribution can be obtained 
by adding the rigid displacement of the unperturbed 
densities, i. e., 

6 Pl (r,qa) = ^>°(r- [RS+<(q<7)]) -p° (r-R»)}. 

a, a 

. (28) 

p® a is the density of the unperturbed ion from Eq. (TTJ 
and the displacement of an ion in Eq. fTS]). 
The results shown in Fig. [7] demonstrate that the non- 
local EPI leads to attractive ( ) and repulsive ( — 

••) regions for electrons mainly at the Bi ions. For 
symmetry reasons some charge is also transferred to 
the silent oxygens in the Of and Of mode, see Fig. 
[7Ji, [TJd and Table |TJ This, however, is not possible 
in Og because all the oxygen ions are vibrating sym- 
metrically towards or away from Bi. Electrons, i. e. 
negative charge, is flowing to the expanded octahedra 
from the compressed ones. Thus, negative charge is 
always accumulated in the larger octahedra and so a 
(dynamical) charge disproportionation of the Bi ions 
is generated. 
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FIG. 7: Contour plot of the nonlocal part of the phonon-induced charge density redistribution from Eq. (|27p for (a) the 
Og mode, (b) the Og 1 mode and (c) the Og mode in model 3. The units are 10~ 4 /a B . Full lines ( — ) describe regions 
of space where electrons are attracted and the lines ( ■ • ) represent regions where electrons are pushed away. 



FIG. 8: Calculated anomalous branches with the oxygen- 
breathing modes Og, Og and Og as end points of the 
A ~ (0, 0, 1), E ~ (1, 1, 0) and A ~ (1, 1, 1) direction, 
respectively. The linetype (• ■ ■ ) denotes the result for the 

RIM (model 1), ( ) additionally includes DF's (model 

2) and the full lines represent the final calculation includ- 
ing DF's and CF's (model 3). 



As mentioned in Ref. |2CJ, if the intersite charge 
transfer in terms of CF's induced nonlocally by the 
OBSM would additionally be spin polarized, as in 
the undoped antiferromagnetic cuprates or the weakly 
doped cuprates with antiferromagnetic correlations of 
sufficient correlation length, corresponding spin fluc- 
tuations are generated (spin-phonon interaction via 
charge transfer). In this way the phonon dynamics 
and the CF's change the spin dynamics and vice versa. 
Such a spin-phonon coupling may be quite strong and 
could lead to a synergetic pairing mechanism in which 
both phonons and spin excitations participate, see 
e.g. Ref. 



26l and references therein. For a qualita- 



tive discussion of a possible interplay between lattice- 
, charge- and spin degrees of freedom in the cuprates 



with the help of the OBSM see also Ref. \M 
The CF's of the Bi6s-, Bi6p- and 02p orbitals <5< K (qcr) 
for 0§, Of and Og are collected in Table Q] Signif- 
icantly larger values are found for model 3 compared 
with the LDA-based model 4. For example in the 0§ 
mode 0.142 electrons are exchanged between the vi- 
brating octahedra in model 3 and only 0.089 electrons 
in model 4. 

For lightly K-doped insulating Ba-Bi-0 besides the 
tilt distortion related to the rotational modes there is 
also a breathing distortion (frozen-in breathing mode) 
of the octahedra. On the other hand, our calculations 
in model 3 for the metallic phase show that strong 
nonlocal EPI of CF-type lead to a large softening but 
a still stable 0§ breathing mode in very good agree- 
ment with experiment. However, the result indicates 
that this type of nonlocal EPI also should play a role 
for the actual breathing distortion in the insulating 
part of the phase diagram. One can interpret the dp- 
patterns in Fig. [7] and the CF's in Table U as an ex- 
pression of a precursor effect for polaron formation. 
Actually, polarons are not expected to form at a high 
enough doping level, because, as mentioned before, 
the system can gain sufficient energy by derealiza- 
tion via the hybridization of the Bi6s and in particu- 
lar the Bi6p states with the 02p hole states created 
by doping. However, such a channel for a lowering of 
the energy by derealization might not be favourable 
in the insulating, lightly doped material due to the 
small number of holes and the tilted structure where 
additionally the 02p to Bi6p hopping is decreased due 
to misalignment of the orbitals and hence the energy 
gain by derealization is reduced. In this case (pola- 
ronic) charge localization by the strong nonlocal EPI 
of CF type might be energetically more favourable and 
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FIG. 9: Same as in Fig. [T]with the calculated results from the model described in the text where additional to the RIM 
dipole fluctuations and charge fluctuations are included. The latter are based on a model where an electronic density of 
states at ef typical for LDA calculations of Ba-Bi-O 2 - is used (model 4). 



responsible for the metal-insulator transition. 
A tendency to polaron formation is also very likely in 
the undoped and weakly doped cuprate-HTSC's, see 
e.g. Refs. I27ll28l . In this context, our modeling of the 
charge response has shown that the delocalized com- 
ponent of the electronic state, i.e., in this case Cu4s, 
is not required to describe the phonon anomalies in 
these " phases" of the cuprates^^^. Polaronic effects 
also occur in the small non-adiabatic region around 
the c-axis of the cuprat-HTSC's where an insulator- 
like charge response and phonon-plasmon mixing is 
founds. 

Our calculations for La-Cu-0 of the orbital resolved 
phonon induced change of the self consistent potential 
8V K felt by an electron, according to Eq. (|21[) . point 
to a hierarchy of energy scales for coupling of the elec- 
trons with the phonons^>2£. There is one energy scale 
related to the OBSM with <5Vcu3d tip to about lOOmeV 
for the Og mode (with a frequency of 79meV) and up 
to about 75meV for the half breathing mode (with 
frequency of 69meV). Both modes have been calcu- 
lated in adiabatic approximation which seems suffi- 
cient for the latter modest. Moreover, we find in La- 
Cu-0 due to weak screening perpendicular to the Cu- 
O plane in a small q space region around the (0, 0, 1) 
direction in particular in the lightly doped materials 
a strong nonlocal, non-adiabatic electron-phonon cou- 
pling of CF type to the charge carriers in the Cu-0 
plane, which are responsible for superconductivity in 
the cuprates, via axial breathing modes of the apex 



oxygen and the La ion. The corresponding values 
of SV K for these modes increase up to about 800meV 
near the phonon-plasmon resonance at low frequencies 
(phonon-like mode, possibly overdamped)^. More- 
over, there is strong coupling at high frequencies be- 
yond the phonon spectrum for the plasmon-like mode. 
Thus, a hierarchy of interaction scales at low energy as 
well as high energy is provided by these nonlocal elec- 
tron phonon coupling effects caused by weak screen- 
ing. Quite recently besides the well known low energy 
"kink" scale (40-80meV), see e.g. Refs. MM, also 
a high energy anomaly in the electron spectrum of 
HTSC's (300-800meV) has been observed by angle- 
resolved photoemission spectroscopy^. See also Refs. 
33.34 for the observation of the high energy anoma- 
lies and an anomalous enhancement of the width of 
the LDA-based CuCVband extending to energies of 
about 2eVM 

As a matter of fact, the electronic structure of the 
lightly K-doped insulating phase of Ba-Bi-0 is not re- 
ally understood. However, as discussed e.g. in Refs. 



35 36 and from our discussion formation of small po- 



laron or bipolarons due to strong nonlocal EPI may 
be an important mechanism for insulating behaviour 
persisting over a wide doping range. 
The mode splitting experimentally found in recent ex- 
periments for the anomalous branches along E and 
A for metallic Ba-Bi-O^ (see the full symbols in the 
figures of the phonon dispersion) cannot be under- 
stood in calculations assuming an ideal cubic per- 
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ovskite symmetry. As discussed in Ref. @ the split 
dispersion suggests some charge inhomogeneity due 
to distortions from the cubic perovskite structure on 
a time scale larger than the time scale of the phonons. 
For example a structural tilt distortion around the cu- 
bic (0, 0, 1) direction recently reported for supercon- 
ducting Ba-Bi-O^l should lead to an increase of the 
OBSM at the M- and i?-point but not at the JT-point, 
as seen in the experiments, because of the decreased 
02p to Bi6p hopping due to misalignement of the or- 
bitals and a corresponding reduction of the softening 
effect via Bi6p. 

Anomalous phonon softening and mode splitting of 
the oxygen bond-stretching phonon along (1, 0, 0), 
most likely due to a (dynamic) charge inhomogeneity, 
very recently also has been observed in La-Ba-Cu-O^. 
As discussed in Ref. a charge inhomogeneity with 
lower symmetry as the " average" translational invari- 
ant charge inhomogeneity considered in our model- 
ing of the electronic state of the underdoped cuprate- 
HTSC's&22 should lead locally to a variation of the 
strength of the nonlocal EPI of CF type and a related 
variation of the phonon frequency of the OBSM. The 
nonlocal EPI may become locally very strong in case 
of a metallic, compressible charge inhomogeneity. 
In our modelling of the electronic state of the under- 
doped cuprates we have incompressible, insulator-like 
regions related to the localized Cu3d states interpene- 
trating a metallic, compressible structure provided by 
the more delocalized 02p states at the O xy sublat- 
tices, where the holes accumulate. Such a correlated 
percolated metallic structure can be thought essen- 
tially to be due to the interplay of the short range 
repulsion (Ud) of the Cu3d states and the long range 
Coulomb interaction together with a gain in kinetic 
energy of the holes by derealization via the 02p 
states. For higher doping in the optimal and over- 
doped state the delocalized Cu4s component of the 
electronic state becomes important in our modelling 
leading to the experimentally observed strong increase 
of phonon softening of the OBSM in the cuprates, see 
e.g. Ref|6J. Physically this means that the perco- 
lated structure melts with increased doping because 
of the enhanced quantum mechanic tunneling via the 
delocalized Cu4s and 02p component in the electronic 
state at the the Fermi energy. Ultimately, in the over- 
doped "phase" where in our modelling the Cu4s com- 
ponent is further enhanced the system should converge 
more and more towards a Fermi liquid. 
The OBSM are a sensible probe for charge inhomo- 
geneities in form of regions with a different doping 
level''. How the frequency would change depends in 
particular on the electronic state of the inhomogene- 
ity which could be a mixture of compressible, metallic 



(higher doping), incompressible, insulating (no dop- 
ing) or partially incompressible (lower doping) regions 
in the crystal. The incompressible and partially in- 
compressible regions (with OBSMs of higher frequen- 
cies) most likely are stabilzed by an interplay of strong 
electronic correlation and strong nonlocal EPI of CF 
type triggering polaronic behaviour, while in the com- 
pressible regions (with lower frequencies of the OB- 
SMs) derealization of the electronic structure changes 
the nature of the electronic state as discussed in Refs. 

H. 

In a phonon-mediated pairing mechanism such a lo- 
cal variation of the nonlocal EPI by inhomogeneity 
is expected to induce simultaneously a corresponding 
change of the pairing strength (gap value) and a vari- 
ation of the frequency of the pairing related modes. 
Usually, stronger renormalized modes (lower frequen- 
cies), i.e., a stronger nonlocal EPI, correlates to a 
larger pairing strength (higher gap value), and vice 
versa. Quite recently such a correlation has been ob- 
served locally in scanning tunneling experiments in 
Bi2Sr2CaCu20g^ and as one possibility attributed to 
electronic-lattice interactions releated to a competing 
electronic ordered state (see, e.g. Ref. |38[ ). 
In the last topic of this paper we present calculations 
for the phonon density of states within model 3 in- 
cluding CF's and DF's, see Fig. [TO] Moreover, in 
Fig. [11] the calculated atom-resolved phonon density 
of states is shown. Concerning the oxygen ion we have 
further discriminated between displacements parallel 
(0||) and perpendicular (OjJ with respect to the Bi-0 
bond in order to demonstrate the strongly anisotropic 
behaviour of the oxygen vibrations in Ba-Bi-O. 
From the figures we see that the frequency range up 
to about 5 THz can be attributed to Ba and Bi vi- 
brations with some contribution from Oj^ including 
the unstable rotational modes. The Oj_ vibrations 
dominate the spectral range from around 5 THz to 
about 11 THz. Beyond this frequency range the high- 
frequency part from 10-15 THz is exclusively governed 
by the Oy vibrations including the von Hove points 
of the OBSM: Of « 13.6 THz, Of w 11.9 THz and 
0§ pa 10.2 THz. The anisotropy between the 0|| and 
O^ vibrations will for example be reflected in the 
Debye- Waller factors. The latter are expected to be 
considerably larger perpendicular to the Bi-0 bond 
in agreement with experiments^ and consistent with 
a tendency towards a structural phase transition via 
rotational modes. Finally, the strong almost singular 
peak at about 15 THz results from the characteristic 
nearly dispersionless branches in Fig. [4] terminating 
at the r point in the TO mode with the highest fre- 
quency. 
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FIG. 11: Partial phonon density of states calulated with model 3. Ba-DOS (a); Bi-DOS (b); ± -DOS (c); Oy-DOS (d). 
The symbol || or _!_, respectively, means parallel or perpendicular to the Bi-O bond. 
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IV. SUMMARY AND DISCUSSION 

Within our microscopic modeling of the electronic 
charge response in linear response theory we have 
studied complete phonon dispersion curves, phonon- 
induced charge density redistributions and charge 
fluctuations, generated by nonlocal EPI, and the to- 
tal and atom-resolved phonon density of states in cu- 
bic metallic Ba-Bi-O. Our calculated results of the 
phonon dispersion curves are in a good overall agree- 
ment with the experiments, in particular the strong 
phonon softening of the anomalous oxygen-breathing 
modes is well described. This was not possible so far 
neither with empirical nor with ab initio calculations 
based on the LDA using a virtual crystal approxima- 
tion. The latter seem to underestimate the kinetic 
single particle part of the density response. Increas- 
ing the electronic PDOS for the Bi6s and Bi6p orbitals 
at the Fermi level as compared with typical values ob- 
tained in LDA calculations has been shown to provide 
the strong softening of the OBSM. This is achieved by 
a characteristic type of local screening on microscopic 
scale in terms of ionic CF's being quite generally at 
work in "ionic" metals with a significant component 
of ionic binding like Ba-Bi-0 or the the cuprate based 
HTSC's. More explicitly, nonlocal coupling of the dis- 
placed oxygen ions in the OBSM essentially to Bi6s- 
and Bi6p CF's, which are favoured by the correspond- 
ing large PDOS of these states at ep, are found to be 
responsible for the softening. 

The effect of the strong nonlocal EPI is also visible 
from the calculation of the nonlocal part of the charge 
density redistribution and the CF's themselves during 
the oxygen vibrations in the octahedra. Accordingly, 
a considerable (dynamical) charge disproportionation 
has been calculated. The patterns of charge rear- 



rangement can be looked upon as a visualization of a 
precursor effect for polaron formation. Small polaron 
and/or bipolaron formation is very likely an impor- 
tant mechanism for the insulating behaviour in lightly 
K-doped Ba-Bi-O. 

The strong effect for softening by including the Bi6p 
orbitals into the electronic state at the Fermi energy 
points towards a change of this state with K-doping. 
This change is accompanied by a corresponding varia- 
tion of the topology of the Fermi surface and possible 
changes of the electron-phonon matrix elements. Both 
effects are important for the electron-phonon spec- 
tral function a 2 F{u>) and the T c value in a phonon- 
mediated mechanism for superconductivity. Such a 
mechanism is supported by the strong nonlocal EPI 
effects found in our calculations. Compared with the 
cuprate based HTSC spin degrees of freedom seem to 
play no role for superconductivity in Ba-Bi-O but the 
strong nonlocal mode-selective coupling effects of lat- 
tice and charge degrees of freedom in terms of ionic 
CF's are present in both, Ba-Bi-O and the cuprate- 
HTSC's. 

From our calculations of the PDOS of the phonons 
we extract that the oxygen vibrations and conse- 
quently also the mean square displacements are very 
anisotropic. There are low-frequency Oj_ vibrations 
perpendicular to the Bi-0 bond which are well sepa- 
rated from the high-frequency 0|| vibrations parallel 
to the bond. 

Quite generally, our investigations for metallic Ba-Bi- 
O and the cuprate-HTSC's point to an interrelation 
between doping induced electronic structure changes 
at the Fermi energy, selective phonon softening via 
strong nonlocal EPI of CF' type and high-temperature 
superconductivity in "ionic" metals. 
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